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Abstract 

Evidence of the absence of the clear electron spin-resonance signal from Ge dangling bonds (DBs) 
at Ge/Ge02 interfaces is explored by means of first-principles electronic-structure calculations. 
Comparing the electronic structures of the DBs at Si/Si02 and Ge/Ge02 interfaces, we found that 
the electronic structure of the Ge-DB is markedly different from that of the Si-DB; the Ge-DB 
states does not position in the energy band gap of the Ge/Ge02 interface while the Si-DB states 
clearly appears. In addition, the charge density distribution of the Ge-DB state spreads more 
widely than that of the Si-DB state. These features are explained by considering the metallic 
properties of the bonding network of the Ge/Ge02 interface and the structural deformation of the 
Ge bulk at the Ge/Ge02 interface due to the lattice-constant mismatch. 
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I. INTRODUCTION 



The length scale of Si-based electronic devices has been continuously decreasing and is 
approaching its technological and physical limits. Thus, a considerable number of studies 
have been conducted to find alternative materials to further improve the performance of 
metal-oxide-semiconductor field-effect transistors (MOSFETs). Among the possible candi- 
date materials, Ge has attracted considerable attention owing to Ge-MOSFETs exhibiting 
higher carrier mobility than Si-MOSFETs. One of the most important issues in the appli- 
cation of Ge-MOSFETs is the accurate control of the properties of the Ge/Ge02 interface 
because it exists even at a Ge/high-A; oxide interface. Electron spin-resonance (ESR) inves- 
tigations have revealed that interfacial dangling bond (DB) defects at Ge/Ge02 interface 
play a different role from those at Si/Si02 interfaces;- a measurable density of DBs of the 
semiconductor surface atoms is not found at Ge/Ge02 interfaces although Si-DBs are ESR- 
active. 

Very recently, Matsubara et al 2 and Hosoi et al— demonstrated that the interface de- 
fect density at Ge/Ge02 interfaces is lower than that at Si/Si02 interfaces when both the 
interfaces are fabricated with no hydrogen passivation treatment. Houssa et al— examined 
by first-principles calculation the density of DBs at the Ge/Ge02 interface as a function of 
oxidation temperature by combining viscoelastic data for Ge02 and the modified Maxwell's 
model, and they addressed the lower defect density at Ge/Ge0 2 interfaces than at Si/Si0 2 
interfaces. In addition, the present authors performed a first-principles calculation of the 
formation energy of defects at the Ge/Ge02 interface during oxidation and found that 
hardly any defects are generated at the interface compared with the number generated 
at the Si/Si02 interface. 5 These theoretical calculations are in agreement with the above 
experiments,-^and the low defect density at the interface may account for the absence of 
DBs observed in ESR experiments, consistent with the conclusions drawn in these studies. 

On the other hand, from the viewpoint of the electronic structure of DBs, Weber et 
a/.- investigated the electronic structure of DBs in Ge bulk by first-principles calculations 
and found that Ge-DBs have an energy level below the Ge valence band edge; Ge-DBs are 
always negatively charged and fully occupied. Broqvist et al- also examined the charged 
states of the DBs in Ge bulk and found two charge-transition levels in the lower part of 
the band gap, indicating the occurrence of three charged states in analogy with the DBs in 
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Si, using the hybrid functional, although these two levels are below the valence band edge 
in the case of the local density approximation. However, Weber et al- have also proved 
that the underestimation of the energy band gap due to the local density approximation 
does not cause the negatively charged state of DBs by employing the Go Wo approximation. 
Houssa et al. 4 reported that the defect state associated with the DBs lies approximately in 
the middle of the Ge energy band gap using a model in which GeC>2 is piled on a Ge(OOl) 
substrate. However, they employed pseudopotentials including a repulsive potential so that 
the energy band gap of Ge agrees with the experimental value, and it is unclear the extent to 
which the deformation of the pseudopotentials affects the atomic and electronic structures 
at the Ge/Ge02 interface. Although it is difficult to adhere to one of these views solely 
owing to the use of approximations in these theoretical calculations, the examination of the 
difference between the electronic structures of the DBs at Si/Si02 and Ge/Ge02 interfaces 
is of importance for clarifying the electronic structure of DBs and for the realization of 
Ge-based devices. 

In this study, we investigate the electronic structures of DBs at the (001)Ge/GeO2 in- 
terface compared with those at the (001)Si/SiO2 interface by first-principles calculations to 
determine the reason for the absence of the distinct ESR signal from DBs at the Ge/Ge02 
interface. Our findings are as follows. The electronic structure of the Ge-DB is notably 
different from that of the Si-DB. The Ge-DB state is occupied at any position of the Fermi 
level, whereas the Si-DB is paramagnetic when the Fermi level is near the midgap. In ad- 
dition, the Ge-DB state couples with the interface state and its charge density distribution 
extends more widely than that of the Si-DB state. 

II. COMPUTATIONAL METHODS 

Our first-principles calculation method is based on the real-space finite-difference approach,- - — 
which enables us to determine a self-consistent electronic ground state with a high degree of 
accuracy using a timesaving double-grid technique.- 1 ^ Intensive studies by various methods 
have shown that the crystalline phase of Si02 can be observed down to ~ 10 A from the 
interface in some cases.— r — Although no crystalline phase of Ge02 has yet been exper- 
imentally observed, we adopt the same atomic configuration for the Ge/Ge02 interface 
to compare the electronic structure of DBs with that of DBs at the Si/Si02 interface. 
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Our interface structure corresponds to that proposed by Kageshima and Shiraishi for the 
(001)Si/SiO2 interface during oxidation.— Figure [1] shows the computational model, con- 
sisting of an a-quartz crystalline oxide layer on a (001) substrate of ~23 A thickness. Since 
Yamasaki et al. 16 claimed that an effective band-gap change from Si to SiC>2 occurs on 
the oxide side within 5 A from the Si/Si02 interface, the thickness of the oxide layer is 
chosen to be ~ 8 A. The size of the supercell is taken to be 2aoX 2aoX 5ao, where ao is 
the experimental lattice constant of Si (5.43 A) or Ge (5.65 A). The periodic boundary 
condition is imposed in the x and y directions, while the isolated boundary condition is 
employed in the z direction. The DBs at the interface are generated using the procedure 
employed by Houssa et al.f a bridging O atom is removed and the lower DB is passivated 
with a H atom. The norm- conserving pseudopotentials^ of Troullier and Martins^ are 
used to describe the electron-ion interaction and are transformed into the computationally 
efficient Kleinman-Bylander separable form,— using the s and p components as nonlocal 
components for O and Si and the s, p, and d components as nonlocal components for Ge. 
Exchange correlation effects are treated by a local spin density approximation, 20 and the 
freedom of spin is considered. The cutoff energy is set to 112 Ry, which corresponds to a 
grid spacing h of 0.30 a.u., and a denser grid spacing of h/3 is taken in the vicinity of the 
nuclei with the augmentation of double-grid points.- 1 ^ Eight-/;; points in the lxl lateral 
unit cell are used for Brillouin zone sampling. Structural optimization is implemented for 
all atoms except the Si (Ge) atoms in the bottom layer and the H atoms on the Si (Ge) 
surface side, reaching a tolerance of the force of F max < 0.1 eV/A. 

III. RESULTS AND DISCUSSION 

Figure [2] shows the densities of states (DOS) of the interfaces without the DB, which 
are plotted by integrating them along the plane parallel to the interface, p(z, E) = 
J \ijj(r, E)\ 2 dru, where r = (x,y,z), tps are the wave functions, and E is the energy of 
the states. For comparison, the DOS of the interfaces without any DBs is also depicted. 
Since we employ a model with a finite-thickness Si (Ge) substrate, the energy band gap 
is overestimated owing to the confinement effect; the energy band gaps of Ge/Ge02 and 
Si/Si0 2 are ~0.20 eV and ~0.62 eV, respectively. However, Gaussian broadening of 0.22 
eV is used to plot Fig. [2] and the band gap of Ge is smeared out by the broadening in Fig. |2j 
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In both the Si/Si02 and Ge/GeC>2 interfaces, the DOS deep inside the substrates is not 
affected by the presence of the DBs. In addition, the peaks of the DOS are observed at 
Ep - 0.3 eV, Ep - 1.0 eV, E F - 1.5 eV, and E F - 2.0 eV, where E F is the Fermi level. The 
states at Ep — 1.0 eV indicated by the black down arrows are related to the Ge- (Si-) O 
bonds since their charge density distributions accumulate between the Ge- (Si-) and O-atom 
layers, and the bunches at Ep — 0.3 eV and accumulating just below the interfaces are the 
contribution of the interface states. When a DB is introduced in the Si/Si02 interface, the 
DOS at Ep — 1.0 eV decreases. Instead, a state appears in the energy band gap, indicated 
by the white down arrow, and the Fermi level is positioned at the state. The midgap state 
is attributed to the n electron of the Si atom because its charge density distribution exhibits 
bunches above and below the Si atom. On the other hand, the peak of the DOS at Ep — 1.0 
eV related to the Ge-0 bonds becomes shallow and the DOS at Ep — 0.3 eV, denoted by the 
white down arrow, increases in the case of the Ge/Ge02 interface. In addition, the Fermi 
level is at the valence band edge of the Ge substrate. This result implies that the Ge-DB 
state of the Ge/Ge02 interface does not have unpaired electrons when the Fermi level is in 
the energy band gap of the Ge substrate, which agrees well with the conclusion of another 
first-principles calculation indicated that the Ge-DB state in Ge bulk is below the valence 
band edge of Ge bulk.- The charge neutrality of the Ge substrate might affect the Ge-DB 
state because one electron is removed from the Ge substrate by the Ge-DB. However, no 
significant differences were found in the DOS when additional Ge layers are attached and 
when one electron is added to the present system. 

We illustrate in Fig. [3] the charge density distributions of the state at the Fermi level of 
the Si/Si02 interface and the state below 0.3 eV the Fermi level of the Ge/Ge02 interface, 
which correspond to the DB states corresponding to the DOS depicted in Fig. [2j In the case 
of the Si/Si02 interface, the charge density distribution of the center— is in agreement 
with the results of a previous first-principles calculation.— The charge density distribution 
of the Ge-DB state is spatially more extended in the lateral direction than that of the Si-DB 
state, and we cannot observe the strong contribution of the ir electron of the Ge atom, which 
can be found in the Si-DB state. Thus, the Ge-DB state couples with the interface states 
of the Ge substrate. Ge is a group IV element between Si and Sn, the latter being metallic, 
and the rutile phase is energetically more stable than the other possible phases of Ge0 2 and 
Sn0 2 23 whereas Si0 2 preferentially adopts the quartz structure, maintaining its sp 3 bonding 
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networks. In addition, the Ge-0 bonds at the interface do not have a strongly preferential 
bonding direction." Thus, the bonding network of Ge/GeC>2 is easily deformed from the sp3 
structure and exhibits more metallic properties than that of Si/SiO^. These characteristics 
of Ge contribute to the smaller energy variation of the Ge-0 bond state after the breaking 
of the bond than that of the Si-0 bond state and the coupling of the Ge-DB state to the 
Ge/Ge02 interface states. These results are in agreement with the experimental result by 
the ESR microscopy- and the absence of the clear ESR signal from Ge-DBs is attributed to 
the significant difference between the electronic structures of Si- and Ge-DBs. Our results 
also indicate that a DB near the interface generates a fixed negative charge, creating serious 
problems for devices that rely on the formation of an electron channel. 

There remains one concern that our main conclusion might be affected by the compu- 
tational models naively. For example, only a limited number of cases have been examined 
in the present work and the Ge-DB states will appear between the valence and conduction 
band edges in some interface atomic structures. Actually, it is difficult to deny the presence 
of the ESR-active Ge-DBs only by the present result. However, the charge density distribu- 
tion of the Ge-DB state is not strongly concentrated around the Ge-DB as shown in Fig. [3] 
and the ESR signal from Ge-DBs is expected to be shallow. 

For the further interpretation of the electronic-structure of the DB states, let us discuss 
about the relationship between the DB states and the deformation of the atomic structure 
at the interface. The interface reconstruction is caused by the lattice-constant mismatch 
between Ge (Si) and Ge02 (SiO"2). The interface Ge atom with the DB is raised up while 
the Si atom is pulled down so that the ir electron forms a P& center in the Si substrate. The 
average variations of the three bond angles formed by the interface atom and two atoms at 
the second interface layer from the ideal diamond structure (109.5 degree) is +5.4 and —8.7 
degrees in the case of the Si/Si0 2 and Ge/Ge0 2 interfaces, respectively. The large variation 
of the bond angles of the Ge/Ge0 2 interfaces is interpreted by the structural properties of 
Ge02 and Si02 under the pressure: O-Ge-0 angle predominantly varies when the Ge02 
is compressed while Si-O-Si angle is deformed in the case of Si02-— Moreover, Haneman 
reported in the study concerning the reconstruction of the Si surface that the sp^ DB bonds 
is stabilized due to the strong contribution of the s electron when the surface atom is raised 
up.— This result also supports our conclusion for the stabilized Ge-DB state and the small 
contribution of the 7r electrons to the Ge-DB state as shown in Fig. [3j 
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IV. SUMMARY 



We have carried out first-principles calculations to examine the electronic structures of the 
DBs at (001)Si/SiO2 and (001)Ge/GeO2 interfaces. We found that the electronic structures 
of the Si- and Ge-DBs are markedly different; the Si-DB state appears at the energy band 
gap of the Si substrate and the position of the Fermi level corresponds to the DB state, 
whereas the Ge-DB state is fully occupied. Moreover, the charge density of the Ge-DB 
state is more widely distributed than that of the Si-DB state. These characteristics are 
explained by the metallic properties of the bonding network of the Ge/Ge02 interface and 
the structural deformation of the Ge bulk at the Ge/Ge02 interface. Our results are in 
good agreement with the finding in ESR measurement at the Ge/Ge02 interface.- Since we 
studied only limited models of interfaces, not all the details of the actual process are included. 
Indeed, our results do not conteract the presence of the ESR- active Ge-DB at the Ge/Ge02 
interface, which have been reported by means of electrically detected magnetic resonance 
spectroscopy.— These findings further motivate the need for a theoretical knowledge for the 
relationship between the atomic structural property of Ge/Ge0 2 intarface and the ESR 
activity of the Ge-DBs. 
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FIGURES 
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FIG. 1. Schematic image of computational model. Black, gray, and white circles represent Si (Ge), 
O, and H atoms, respectively. The rectangle enclosed by broken lines represents the supercell. 
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FIG. 2. Distributions of DOS integrated on plane parallel to the interface as functions of relative 
energy from the Fermi level, (a) Si/SiC>2 without DBs, (b) Si/Si02 with DBs, (c) Ge/Ge02 
without DBs and (d) Ge/Ge02 with DBs. Zero energy is chosen to be the Fermi level. Each 
contour represents twice or half the density of the adjacent contour lines, and the lowest contour 
is 1.45 xlO -4 e/eV/A. The dashed and dotted lines represent the vertical positions of Si (Ge) 
and O atomic layers, respectively. The right arrows denote the position of the Si and Ge atoms 
forming DBs. Gaussian broadening of 0.22 eV is used and the band gap of Ge is smeared out by 
the broadening. 
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FIG. 3. Contour plots of charge density distributions of DB states. The planes shown are along 
the cross section in the (110) plane including DBs. (a) Si-DB state at Ep and (b) Ge-DB state at 
Ep — 0.3 eV. The lowest contours are 6.75 x 10 -3 e/A 3 and the subsequent contour lines represent 
values larger than 6.75 x 10~ 3 e/A 3 . The meanings of the symbols are the same as those in Fig. [TJ 
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